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The partition function of the two-dimensional Ising model with zero magnetic field on a square 
lattice with m x n sites wrapped on a torus is computed within the transfer matrix formalism in an 
explicit step-by-step approach inspired by Kaufman's work. However, working with two commuting 
representations of the complex rotation group SO(2n,C) helps us avoid a number of unnecessary 
complications. We find all eigenvalues of the transfer matrix and therefore the partition function in 
a straightforward way. 

PACS numbers: 05.50.-|-q, 68.35.Rh, 05.70.Ce, 64.60.Cn 



I. INTRODUCTION 

Since Onsager's solution [1] in the transfer matrix ap- 
proach [2] of the two-dimensional Ising model [3] with 
vanishing magnetic field and its subsequent simplifica- 
tion by Kaufman [4] , there have been a number of related 
as well as alternative solutions, see e.g. Baxter's book on 
exactly solved models in statistical mechanics [5] and ref- 
erences therein. Among the transfer matrix solutions are 
the ones by Schultz, Mattis and Lieb, by Thompson, by 
Baxter and by Stephen and Mittag [6]. Nevertheless, the 
author of this work feels that there is still room for a nice 
and straightforward solution. A much abbreviated ver- 
sion of this work, sketching only the new aspects may be 
found in [7], while here we give a self-contained account 
of our approach with all relevant details. 

We study the two-dimensional Ising model with zero 
magnetic field on a square lattice with m rows and n 
columns subject to toroidal boundary conditions. The 
transfer matrix is expressed in terms of the generators of 
two commuting representations of the complex rotation 
group S0(2n, C). These representations naturally arise 
from projected bilinears of 2" x 2" spin matrices. Con- 
servatively speaking, we reduce Kaufman's approach to 
its essential steps, avoiding in particular the doubling of 
the number of eigenvalues of the transfer matrix and sub- 
sequent rather involved arguments for the choice of the 
correct ones. We also need not investigate the transfor- 
mation properties of the spin matrices. In our opinion, 
the value of this work is not that it gives the most concise 
solution, but that it provides a clear line of arguments 
and requires only rather basic tools from mathematics 
and mathematical physics. 

We proceed in a step-by-step fashion, exhibiting also 
those steps that closely parallel the corresponding steps 
in Kaufman's treatment [4] or Huang's corresponding 



textbook write up [8] . On the one hand this is necessary 
because of our slightly different (but, in our opinion, more 
natural) conventions, in particular for 70. On the other 
hand, we would like to keep this work as self-contained 
as possible. 

Finally, this work may help sharpen the view on the 
difficulties encountered when attempting to solve the 
two-dimensional model with magnetic field or the three- 
dimensional model and might be useful for simplifying 
other problems in statistical mechanics for which an ex- 
act solution is attempted. 

Our notation is in the spirit of [8] , with sans serif cap- 
itals reserved for 2" x 2" matrices. The structure of this 
work is as follows: In section II the model and its trans- 
fer matrix T are defined. In section III we express the 
transfer matrix in terms of 2" x 2" spin matrices Xj^, 
Yi^, Z^. A rescaled transfer matrix V is defined whose 
eigenvalues are, up to a trivial factor, the eigenvalues of 
T. In section IV, we define further spin matrices Vi, and 
two commuting projection classes J^^ and of their 
bilinears. After investigating the relevant properties of 
the J^^, we express V in terms of them. In section V, 
we introduce 2n x 2n matrices Jap whose algebra is iden- 
tical to that of J^fl and define matrices in terms of 
the Jaf3 such that the relation between and Jajs is 
closely related to that between V and J^^. The are 
subsequently diagonalized. In section VI, the analogy 
between and V is exploited for the diagonalization of 
V. The eigenvalues of V and the partition function are 
found explicitly. In section VII we comment on the tech- 
nical difficulties encountered when attempting to solve 
the two-dimensional model with magnetic field or the 
three-dimensional model. 



II. THE MODEL AND ITS TRANSFER 
MATRIX T 
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We work with a square lattice with m rows and n 
columns and consequently m x n sites. The energy is 
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FIG. 1: Lattice of the two-dimensional Ising model. The 
parameters a and h correspond to the interactions in the re- 
spective direction indicated. Rows 1 and m+1 and columns 1 
and n+1 are identified, respectively, i.e. the model is wrapped 
on a torus. 



where tt^ = {s^i, . . . , s^„} for /i = l,...,m and T is 
defined by its elements (s„+i = si), 

n 

(7r|T|7r') = JJ exp(as^s'^ -|- hSvSv+i). (6) 

We can split T into a product of two matrices 

T = V,V;, (7) 
defining and V(, by their elements 



L/=l 



(8) 



and 



given by 



E = Ea + Eb, 



(1) 



(7r|V6|7r') = Y[ Ss^si exp{bSt,Sv+i). 



v=l 



(9) 



-l3Ea = a'^ SiSj = a'^^Sf^^Sn+i^^, (2) 

{ij)a ,i=li/=l 



-(3Eb = b^SiSj = b^^Sfj,^Si^^^+i, (3) 

{ij)b /i=li/=l 

with a = —PJa, b = —(3Jb and f3~^ = kT, where Jo 
and Jft are temperature-independent interaction energy 
parameters which by convention are negative (other sign 
combinations do not change the partition function, as 
can be seen by appropriate redefinitions of half of the 
signs of the spin variables). Consequently, a and b are 
real and positive. {ij)a and {ij)h indicate summation 
over neighboring spins in the respective directions. Wc 
identify rows 1 and m+1 and columns 1 and n+1, i.e. the 
lattice is wrapped on a torus, see also Fig. 1. The s^i, 
can take the values ±1. The partition function is then 
given by 

Z(a,6) = ^---^exp(-/3E) 



= I^'-'I^exp 



.fi=l jy=l 



= • • • XI n n *5^p(«sm!^'Sm+i,>' + bSfi^Sfi^^+i). (4) 

Sll Smn JU=1 y=l 

It can be expressed with the help of a 2" x 2" transfer 
matrix T, 

Z{a, b) 

= ^{'^l\'^Wm){'^m\'^\'^m-l) ' ' ' (7r3|T|7r2)(7r2|T|7ri) 



III. T AND SPIN MATRICES X^,, Y^, 



With the help of the Pauli matrices 



Cx = 



1 



CTy = 



0-z 



-1 

(10) 

and the 2x2 unit matrix 1, define Hermitian 2" x 2" 
spin matrices by the direct products 

X^=l®---®l®CTx®l®---®l, (11) 

^ V ' ^ V ' 

v—1 n—u 

and analogously for and Z^. For any v and v' ^ they 
obey the commutation relations 

[X.,X,,] = [Y,,Y,.] = [Z,,Z,,] = 0, (12) 
while only for v ^ v' holds 

[X., Y.'] = [Y,, Z,.] = [Z,, X,.] = 0. (13) 
For any v we have 

(14) 



= Y^ = Z^ = 1, 



(5) 



{X,, Y4 = {Y,, Z4 = {Z,, X4 = 0, (15) 

X,Y, =iZ„ Y,Z, = iX„ Z,X, =iY„ (16) 

where 1 is now the 2" x 2" unit matrix. 
Define a by 

sinh(2a)sinh(2a) = 1, (17) 

so that a > 0, 

tanh o = exp(— 2a), tanha = exp(— 2a) (18) 
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and 



e " e ' 
Then we can write 

V; = [2sinh(2a)]"/2V„ 

with 



„ l+a I = [2sinh(2a)]V2exp(aa,). (19) 



(20) 



and 



Va = ]^ exp (aXj,) , 



Vfc = exp (6Zi.Zi.+i) , 



(21) 



(22) 



where we have identified Z„+i = Zi. The transfer matrix 
(7) may then be expressed as 



T= [2sinh(2a)]"/2V6V„. 



(23) 



Due to the c;ychc property of the trace, we may rewrite 
the partition function (4) as 

Z{a, b) = [2 sinh(2a)]'""/2TrV™, (24) 

where V is defined by the Hermitian matrix 

V = V„/2V6V„/2 (25) 

with 



V„/2 = n exp(aX,/2), 



(26) 



i/=i 



so that V^^2 = Va ■ If Afe are the 2" eigenvalues of V, we 
have 



2" 



Z{a, b) = [2 sinh(2a)]'""/2 ^ a^. (27) 



fe=i 



Our task is therefore to find the eigenvalues of V. 



IV. SPIN MATRICES AND ALGEBRAS OF 
THEIR PROJECTED BILINEARS J^^ 



Define the 2n matrices {u = 1, . . . ,n) 

^2v-l = Xi • • • Xj,_iZy, 

^21^ = Xi • • • Xj,_iYy, 



(28) 
(29) 



which obey 

{r^,r,} = 2d^^. (30) 

Define further the matrix 

Ux = Xi---x„ = i"rir2---r2„, u2 = i, (31) 



which anticommutes with every f^, 
{^,Ux} = 0. 



(32) 



Now we can write for the matrices appearing in the ex- 
ponents of (22) and (26) 

X. = -^[r2.,r2.-i], u = l,...,n, (33) 



Z;yZ^_|_l 



;[r2,.+i, r2,.], i/ = l,...,n- 1,(34) 



z„Zi = -Ux[ri,r2„]. 



(35) 



So far our treatment has been rather similar to Huang's 
write up [8] of Kaufman's approach [4] . Our subsequent 
treatment rests on the observation that the formulation 
(25) of V with (22) and (26) involves only the product 
Ux of all and bilinears ^a^js, see (31), (33)-(35). This 
will allow us to find two commuting algebras of projected 
bilinears of the f^. 

Define the projectors 



P±^-(l±Ux), 

for which hold 

P+-|-P- = l, P+P+ = P+, p-p- = P 



(36) 



p+p- = p-p+ = p±y^ ^ UxP=^ = ±P=t ^^''^ 



and 



[P^,r«r^] = o. 



With their help, define the matrices 



Ja/3 — —-^[^a, T/s], J^^ — P^Ja/3, 



(38) 



(39) 



so that 



Since 



(41) 



there are n(2n — 1) such independent matrices of each 
kind J^^ and J~^. It is straightforward to show that 
their algebra decomposes into two commuting parts, 

[^lp,J^s]=0, (42) 
which obey identical algebras 

Next note that with (33)- (35), (39) and (40) we can 
write 

X. = 2(J+^2.-i+J2;..2,.-i). '^ = l,...,n, (44) 
Z,Z,+i = 2(J++i2,+J2-+i2J: 1^= l,...,n-l,(45) 

Z„Zi = -2Ux(J+2n+Jr,2„) = -2(Jt2n--lr,2n)- (46) 
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This allows us to express from (26) and from 

(22) in terms of the J^^, 

n 

Va/2 = n exp[a(J2+.,2.-i + J2".,2.-i)] = Kn^^l^ (^7) 



with 



and 



^a/2 = n "Xp(ajJ^_2i.-l)> 



(48) 



Vb = exp[-26(J+2„ - Jr,2j] 

n-l 

X ]J exp[26(J+ +i_2^ + -'2^+i,2>.)] 



(49) 



with 



V± = exp(T26J±2j n exp(26J±,+i,2J- (50) 



The rescaled transfer matrix V defined in (25) reads then 
V = V+V" (51) 

with 

V^=Vj2V^V^/2' [V+,V-]=0. (52) 



V. THE MATRICES AND THEIR 
DIAGONALIZATION 

In this section we formulate matrices whose defini- 
tions are similar to the expressions of V='= through (48), 
(50) and (52). However, the have convenient peri- 
odicity properties which allow for explicit diagonaliza- 
tion. Our treatment closely parallels that of Kaufman 
[4] as written up by Huang [8]. However, we will be 
more explicit and also carefully analyze some subtleties 
of the analysis concerning the representation of the rele- 
vant similarity transformation matrices S± and the sign 
of 7o, see below. 

Define N x N matrices Ja^ by their elements 



(53) 



where Greek and Latin indices run from 1 to N . Since 

Jap = —Jpa, there arc N{N — l)/2 such independent 
matrices. As can be easily checked, they obey the algebra 

[Ja0, J-iS] = i{Sa-fJf3S + Sf3sJa-f — (^q^ J/37 — bfi^Ja_&), (54) 



which is, if we set N = 2n, identical to the algebras (43) 
of J^^ and J~^. Now consider the matrices 

S = exp{iCal3Jal3), (55) 



for which holds 



S"^ = exp(iCa/3 Jc/s)^ = exp(iCa/3 Jj^) 



(56) 



Since 



(det Sf = det S'det S'^ = det Sdet S'^ = 1, (57) 

and because S is smoothly connected to the unit matrix, 
we have detS* = 1. For real parameters Cap, the S are 
also real and form the group SO(A^) oi N x N orthogo- 
nal matrices with unit determinant. The algebra (54) is 
therefore called the Lie algebra of SO (A''). Here, we let 
the Cq,/3 be arbitrary complex numbers, so the matrices 
S form the group SO(A'^, C) of complex N x N matrices 
with 



=s- 



detS 



Define the S0(2n, C) matrices 



with 



K/2 = JJ exp(aJ2^,2!y-i) 



(58) 



(59) 



(60) 



i/=i 



and 



n-l 



exp(=F26Ji,2n) exp(26J2i.+i,2i/), (61) 



in analogy with V±, V^^2 and Vf in (52), (48) and (50). 



and V± 

Since a and h are real, the matrices V^/2, 



and F± 

are not only orthogonal, but also Hcrmitian, so the 
have only real eigenvalues and in each case a complete 
set of orthonormal eigenvectors. Since the algebras of the 
three sets of matrices J^^, and Jap are identical and 
because will turn out to be explicitly diagonalizable, 
we will be able to find all eigenvalues of V. 
The structure of 14/2 and V^"^ is given by 



a/2 



Ra 
















V 



Ra 



\ 














(62) 



with 



Ra 



cosh a i sinh a 
-isinha cosh a 



(63) 



5 



and 



with 



/ cosh(26) 





i?26 



±?:sinh26\ 







V T«sinh26 



R2b 



(64) 



cosh(26) / 



R2b 



cosh 26 isinh26 
-isinh26 cosh 26 



(65) 



Carrying out the matrix multiphcation in (59) is straight- 
forward and gives 



The equations with through z"~^ on the right hand 
side are identical, which leaves three independent equa- 
tions 

(A + T -2""^-B^)w = Au, (76) 
{A + zB + z-'^B'')u = Xu, (77) 
{Atz^-''B + z-'^B'i)u = Xu. (78) 

These equations can simultaneously be solved when 

z" = Tl. (79) 

since then they become identical, 

{A + zB + z-^B'')u = Xu. (80) 

The sign =F in (79) is associated with V^. Altogether, 
there are 2n values of z that solve (79), 



/ ^ 


B 


• 


• 




^^ = 6"*^/", fc = 0,...,2n-l. 


St 


A 


B ■ 


• 








St 


A ■ 


• 

















(66) Even k lead to 2;" = +1 while odd k lead to 2;" 








• 


• A 


B 


fc = l,3,5,...,2n-l forF+, 







• 


■ St 


A J 


fc = 0,2,4,...,2n-2 for F". 



(82) 



with 



and 



A = cosh 26 



cosh 2a i sinh 2a 
— isinh2a cosh 2a 



B = sinh 26 



— h sinh 2a —i sinh^ a 

2-1 - 

i cosh a — 2 sinh 2a 



(67) 



(68) 



Let us make the following ansatz for the normalized 
eigenvectors of V^, 



( zu \ 

2 



(69) 



where 2; is a complex number and u is a two-component 
vector that we assume to be normalized, u'^u = 1. The 
condition 



V^ip = Xip 
leads to the n equations 



{zA + z'^Btz"B^)u = zXu, 

{z'^A + zJB + zB'^)u = z'^Xu, 
{z^A + z^B + z'^B'')u = z^Xu, 



{z"'-^A + z^'B + z''-^B'')u 
{z''AtzB + z'^-^B'^)u 



z'^-'^Xu, 
z^-Xu. 



(70) 



(71) 
(72) 
(73) 

(74) 
(75) 



For each /c = 0, . . . , 2n — 1, we still have to find the asso- 

\^ and 



dated two eigenvalues AI and Xi and their corresponding 



eigenvectors uj, and uj. in (80), i.e. in 

with the Hermitian matrix 

Mk = A + e''^'^/"^ + e-'''^/"Bt = 



with 



and 



nk 



Explicit evaluation gives 

detMfc = 4 - Okol = 1, 



(84) 



(85) 



dk = cosh 2a cosh 26 — cos — sinh 2a sinh 26 (86) 

n 



Ok = — sin — sinh 26 

n 

-|- i f sinh 2d cosh 26 — cos — cosh 2d sinh 26 ) . (87) 



(88) 



implying [up to an overall sign, fixed by trMfe = A^+A^ = 
2dk > 0, see (86) above] eigenvalues of the form 

AT=e+^^ Ai = e-^'= (89) 
with real 7^. The value of jk can be found from 

trMfe = e+~"' + e^'"' = 2cosh7fe. (90) 
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FIG. 2: for n = 3 and n = 25 as a function of a in the model 

with b — a. The sohd hnos from the top are 7„, 7n-i, • • • , 
70, while the dashed hnes from the bottom are 7^^ T^— i? 

— 7o- Apart from ±70 and ±7,,, all ±7^ are two- fold 
degenerate. For n — + 00, the lines merge into two continua 
above and below 0. 



so that 



nk 



cosh 7fe = cosh 2a cosh 2b — cos — sinh 2a sinh 26, (91) 

n 

and from (88) follows then that we may write 

Ok = ie'^" sinh7fe (92) 

with some phase 5^ to be considered further below. 

If 7fc is a solution, then also — 7^ is one, but this has al- 
ready been taken into account in (89). Let us fix the sign 
of 7fe by defining 7^ = 2a for 6 = and then analytically 
continuing to other values of 6. For k = 1, . . . ,2n — 1, 
this means 7^ > 0. On the other hand, for 70 this means 



7o = 2(a-6). 



(93) 



Our sign convention for the 7^ and in particular for 70 
will allow us to treat all 7fc on an equal footing for both 
a > b and a < b, i.e. irrespective of the temperature. 
It is obvious from (86) that for fc = 1, . . . , n — 1 



(94) 



7fe = l2n-k, 

t I 

i.e. the eigenvalues with k = 1, • . • , 
pairs. Since for < A; < n [extending for the moment 



1 occur in 



(91) to non-integer values of k], 

d-jk TT sinh 2o sinh 26 . nk ^ 
dk nsinh7fe n ' 



we have 



This means 



< |7o| < 71 < ••• <7„. 



(95) 



(96) 



(97) 
(98) 



< +jo < 71 < • • • < 7rn a > 6, 
< -70 < 71 < ■ ■ ■ < 7n) a <b. 

We have plotted two examples with a = 6 in Fig. 2. 
Comparison of (87) with (92) gives 

cos 5k sinh 7/j = sinh 2d cosh 26 — cos — cosh 2d sinh 26, 

n 

(99) 
(100) 



sin 6k sinh7/j = sin — sinh 26. 



The Sk are smooth functions of a and 6. For k = 0, (99) 
and (100) give 



cos 60 sinh 70 = sinh 2 (a — 6) , 
sin 5o sinh 7o — 0, 



so that with (93) we have 

^0 = 0, 



(101) 
(102) 



(103) 



while from (94), (99) and (100) follows for fc = 1, . . . ,n 



'J2n-k 



-5k, 



implying 



^n = 0. 



(104) 



(105) 



The normalized eigenvectors of Mk may now be written 



as 



1 

71 



V2{ 



e 2- 



(106) 

as may be checked by inserting them into (84) . The eigen- 
vectors ip^ of V^, defined according to (69) and (81) by 



^1' = ^ 



2ikTT /n^. Ti. 

A; 



(107) 



are all orthogonal. For different eigenvalues this is clear, 
while for the degenerate eigenvalues 7^ and 72ra-fe, k = 



7 



1 , . . . , n — 1 , this is easily seen by explicitly computing 

and 



the scalar product of ipl^ and ipJ^ 



/ Tit ; Ti 



,Tit„,Ti 



n 



— ^"'k "■2n-k 

= 0. 



1=1 

■j^ g— 27ri(n+l)/c/n 

2 ^ — 2'Kik/n 



1 

(108) 



From here on, we will be considerably more explicit 
than [4] (the following discussion is omitted altogether 
in [8]). This is because, without the following steps, the 
solution is incomplete. 

Let us define 



R 



2n-3' ■ 



and 

(V'd,V'0.V'2>V'2n-2>V'4>V': 



■ >^L-3'V'iV'L-i.V'i) 

(109) 



.V'L-4>V'4>^i'2„-2.V'2)- 

(110) 

The -R_|_ are then unitary and diagonalize according 
to 



R+V+R+^ 



-71 



=+73 



-73 



o-|-72n-l 



-72n-l 



(111) 



and 



-1-70 



0-I-72 



o-|-72n-2 



V 



where (94) has been used. 



,(112) 



2TI-2 J 



The diagonal form above is not the most useful for the 
following considerations. Let us instead apply a further 
similarity transformation to obtain again a form that is 
part of '&0{N,C). First, there is a more useful way of 
writing R^ . For this purpose, define the 2 x 2-matrices 

1 / expi(^ + |) -iexpi(-^-l) 
^V-zexpi(^-l) expi(-^ + |) 



and organize them as 

/ ^11 D:ii 



and 



RZ^ = 



1 



Di2 D32 

\ Din D^n 

D02 D22 



D2n-l,l \ 
D2n-1,2 



D2n-2,1 \ 
D2n-2,2 



(113) 



(114) 



(115) 



\ Don D2n ■ ■ ■ D2n-2,n J 

Now, define Rx and its inverse by 



R 



±1 

X 



/ exp(=F«7rc7x/4) 



exp(=Fi7rcra;/4) 



exp(=Fi7rcra;/4) / 
(116) 



with 



W 1 



exp(Tmc7,/4) = — i^^^^^j (117) 

and apply the further similarity transformation 

= RxR±V^Rg^Rx^ = S±V^S^^ (118) 
to (111) and (112). 



This gives on the one hand 
/ exp(-7{i}C7y) 



exp(-7{3}(7y) 



V 



exp 



exp(-7|2„-i}0-y) / 



-7{3}0-j; 



8 



(119) 



while on the other hand we have 



1 



13 



f Du D 

\ Dni Dn3 

/ Dio D12 



-Dl,2n-1, \ 
D2,2n-1 



D. 



n,2n— 1 



(120) 



S_ ^ — R_ ^ -^x^ 



1 



D20 D 



22 



£'l,2n-2 \ 
-D2,2n-2 



\ DnO Dn2 ■ ■ ■ Dn,2n-2 / 



with 



cos(^ + |) -sin(^ + |) 
sin(^-l) cos(^-l) 



(121) 



(122) 



Let us show now that the S± are elements of S0(2n, C). 
Since R± and Rx are unitary, the S± are also unitary. 
Since the S± are also real, they are orthogonal. We still 
need to show that det5± = 1. Since S± are orthogonal, 
we have 

(det S±f = det S± det 5^ = det S± det = 1 (123) 

and the question is only if dot S± = +1 or dot S± — —1. 
Since S± are analytic in b, we may work with 6 = and 
therefore ^fc = and then analytically continue to non- 
zero b, upon which dot S± cannot change discontinuously 
and therefore not change at all. Wc have 



det5±^ = det R^^S^^Rx = det R^^R^^, 
where for ^fe = 



(124) 



D 



13 



Rx^R+^ 



1 



D21 D23 

\ Dnl Dn3 



Dl,2n-1 ^ 
-D2,2n-1 



Dn,2n-1 / 



(125) 



R^^R_^ 



1 



D20 D22 ■ ■ ■ D[ 



l.ln-l 



(126) 



with 



\ DnO Dn2 ■ ■ ■ -D„,2n-2 / 

exp{+ikl'7r / n) 

exp(— ifc^7r/n) 



(127) 



For the purpose of computing their determinants, we may 
reorganize R'^R^ into block matrices 



Rt 



R4 



(128) 



such that all cxp{+iklTT / n) are in the nx n matrices i?^ 
and all exp(—ikln/n) in the nx n matrices R^, keeping 
the ordering of the exp{+iklTr / n) among themselves and 
of the cx];>{—ikl'K /n) among themselves. Such a reorgani- 
zation involves equal numbers of exchanges of rows and 
columns, respectively, and does therefore not change the 



determinant of R-^^R±.^. Also, we have R^. 
therefore 



and 



detSg^ = det Rx^Rg^ 



det R^ det R^ 



which leaves only 



|deti?i|2 > 0, 



det5± = 1. 



(129) 



(130) 



This means that the matrices S± are elements of 
S0(2n, C) and can be represented as 



S± 



exp(ic^^Ja/3) 



(131) 



with certain unknown complex parameters c^^. 

Let us summarize the crucial result of this section: 
Applying similarity transformations = S±V'^S^^ 
with certain matrices S± = exp(zc^^ J^/s) to V"^ = 
Va/2V^Va/2 with Va/2 = lYl=i exp{aJ2u,2u-i) and = 
exp(=F26Ji,2n) n"=i sxp(2W2,y+i,2iy), we obtain = 
exp (l]"=i7{i-z^}'^2!y,2iy-i) with "fk as defined in (91) 
and the subsequent sign convention. This is what we will 
use in the next section to determine all eigenvalues of V 
and therefore the partition function Z{a,b). 



VI. DIAGONALIZATION OF V 

Now use the same parameters c^^ from the last section 
to define the 2" x 2"-dimensional transformation matrix 

S = S+S_, S± = exp(ic±^J±^), (132) 

and write 

Vs = SVS"^ = S+V+S+^S_V-SI^ = V+Vg. (133) 
The factors 

V| = exp(^c±^J±^)V± ,V±V± 2 exp(-zc±^ J±^) (134) 
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have the same structure as 



Vg = exp(ic±^ J„/3)V;/2V;=^V;/2 exp(-ic^^ J^^j) (135) 

from the last section. Now imagine using the Baker- 
Campbell-Hausdorff formula [9] 



exp(A) exp(S) = exp ( A + S - A] 



+ -{[A,[A,B]] + [B,[B,A]]} + - 



(136) 



to work out all products of exponentials in . Since the 
-'a/3' -'a/3 ^^"^ obcy identical algebras, the result is 



Vf = exp ^7{i-j}J2.,2.-i > 



(137) 



so that 

Vs = exp 

= exp 



.1^=1 
1 " 

-5] 72.-1(1 + Ux)X. 
i/=i 

+ lE ^2.-2(1 -Ux)X. 



(138) 



with the same as defined in (91) and the subsequent 
sign convention. 

To diagonalize Vs, define another similarity transfor- 
mation 

Vy = RyVsRy^ (139) 
with Ry and its inverse given by 

n 

R±i =2-"/2 J|(l±iY,). (140) 
i/=i 



Since 

RyXyRy^ = Zy, 
this transformation takes Vs into 



(141) 



Vy = RyVsRy^ = exp 



1 " 

-5^72.-i(l + Uz)Z. 



1 " 

+ lE ^2.-2(1 -Uz)Z. 



with 



U2 



Zi--Z„. 



(142) 



(143) 



The matrix Vy is diagonal, but we still have to deter- 
mine its elements. Uz is a diagonal matrix with elements 
-|-1 and —1 occurring in equal numbers. For each element 
holds: If an oven (odd) number of Z,^ provides a factor 
— 1, the matrix element of Uz is +1 (—1). This means: (i) 
A matrix element of (1 -|- Uz)/2 is 1 (0) if an even (odd) 
number of Z^ provides a factor —1: (ii) A matrix clement 
of (1 — Uz)/2 is 1 (0) if an odd (even) number of pro- 
vides a factor —1. It follows that the 2" eigenvalues of V 
split into 2"~^ eigenvalues of the form 



exp 



^E(±)72.-ij 



(144) 



and 2" ^ eigenvalues of the form 



(145) 



where in the first (second) case all sign combinations with 
an even (odd) number of minus signs occur. This is re- 
flected by the indices "e" and "o" in our result for the 
partition function. 



i/=i 



Z(a,6) = [2sinh(2a)]'""/2 
= i[2sinh(2a)]'""/2 

C 71 n 

H [2cosh(^72fc_i)] + J] [2sinh( 



-72fe-l 



n n 

] + J] [2 cosh (^^2^-2)] - n [2sinh (y 72^-2) 



(146) 



The last term within the braces has a sign differing from 
that in [4]. This is due to our different sign conven- 
tion for 7o. The eigenvalues of T are of course obtained 
by multiplying (144) and (145) with the trivial factor 



[2sinh(2a)]"/2. 

The results for the eigenvalues of T and the partition 

function are the starting point for the analysis of the 
thermodynamic properties of the two-dimensional Ising 
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model, the most interesting case being the thermody- 
namic hmit m,n ^ oo. Such analyses can now proceed 
as usual (see e.g. [1, 4, 8]) and will not be repeated here. 

VII. REMARKS ABOUT MAGNETIC FIELD 
AND THREE DIMENSIONS 

Let us briefly remark on the difficulties encountered 
when trying to solve the two-dimensional model with 
magnetic field or the three-dimensional model in our ap- 
proach. In the first case, we need besides and V;, 
another matrix 

n 

Ve = n exp(cZ,), (147) 

1/=! 

and in the second case another matrix 

n 

Vd^l[cxp{dZ,Z,+n'), (148) 
v=i 



if the three-dimensional model is obtained from the two- 
dimensional one by letting sites v and ly + n' for v = 
1, . . . ,n interact with each other (assuming we have cho- 
sen n such it can be divided by n'; the lattice is then 
m X n' X n/n' and has in the n' and n/n' directions the 
character of a screw on a torus, while in the m direction 
it remains periodic). For the three-dimensional model 
with magnetic field, we need both Vc and V^. 



Both the Zi, and the Z^Z^+n' (with n' > 1) are not 

part of the algebra of J^^ and J~^. In fact, it turns 
out that amending the Lie algebra with either of these 
classes of matrices makes the number of elements in the 
new algebra grow proportional to 2" instead of n^. This 
eliminates the possibility of using an extended version of 
the algebra of matrices Ja^j with side length growing pro- 
portional to n. But this was essential for the convenient 
periodicity properties of the matrices V^. 
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